Phase coherence of conduction electrons below the Kondo temperature 
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We have measured the phase decoherence rate, r^ 1 of conduction electrons in disordered Ag 
wires implanted with 2 and 10 parts per million Fe impurities, by means of the weak localization 
magnetoresistance. The Kondo temperature of Fe in Ag, Tk ~ 4 K, is in the ideal temperature 
range to study the progressive screening of the Fe spins as the temperature T falls below Tk- The 
contribution to t^ 1 from the Fe impurities is clearly visible over the temperature range 40 mK - 10 K. 
Below Tk , t^ 1 falls rapidly until T/Tk ~ 0.1, in agreement with recent theoretical calculations. 
At lower T, t^ 1 deviates from theory with a flatter T-dependence. We speculate that this latter 
behavior is due to incomplete screening of the s=2 Fe impurities by the conduction electrons. 
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The Kondo problem is a paradigm many-body prob- 
lem in condensed matter physics. In recent years, Kondo 
physics has been observed in semiconducting quantum 
dots and carbon nanotubes, while Kondo lattices play 
a major role in some strongly-correlated materials |2j. 
The original problem that motivated Kondo's 1964 paper 
was the increase in resistivity at low temperature of 
metals containing magnetic impurities. Although Kondo 
solved this puzzle, his perturbation theory diverged in 
the limit of zero temperature. Wilson showed that, 
at temperatures far below the Kondo temperature Tjc, 
the magnetic impurity forms a spin-singlet with the sur- 
rounding conduction electrons and behaves as a non- 
magnetic scatterer with cross section given by the uni- 
tarity limit. 

The brief history given above might leave the impres- 
sion that the the behavior of dilute magnetic impurities 
in metals is completely understood. That is, however, not 
the case. An important aspect of the Kondo problem con- 
cerns the distinction between elastic and inelastic scatter- 
ing. This distinction is extremely important in the con- 
text of quantum transport and mesoscopic physics, where 
it is known since 1979 that elastic scattering from static 
disorder preserves quantum phase coherence, whereas in- 
elastic scattering destroys it. The magnetic impurity 
contribution to the conduction electron phase coherence 
rate, tI 1 , was first measured explicitly by two groups in 
1987 |a,l3 j an d nas received renewed attention recently 
0, 01 in the context of the debate over zero-temperature 
decoherence in disordered metals 0, . Until very re- 
cently 0,0, however, there was no theoretical expres- 
sion for the temperature dependence of the inelastic scat- 
tering rate valid for T not too far below 7V, and very 
little reliable data in that temperature range |13| . 

The most reliable estimates of t^ 1 come from anal- 
ysis of low-field magnetoresistance data in the context 
of weak-localization theory. In disordered metals with- 
out magnetic impurities, t^ 1 is dominated by electron- 
phonon scattering at temperatures above about 1 K, and 
by electron-electron scattering at lower T. In the pres- 



ence of magnetic impurities, contains the additional 
contribution j m , which peaks at T w Tk- In order to 
observe this peak in 7 m while keeping the magnetic im- 
purity concentration low enough to avoid interactions be- 
tween impurities, one must choose a system with Tk be- 
low about 10 K; otherwise t^ 1 is dominated by electron- 
phonon scattering. In order to acquire data far below 
Tk, however, it is important to keep Tk as high as pos- 
sible. The optimal range for Tk is a few Kelvins, which 
is achieved with Fe impurities in Ag |l4| . 

We fabricated Ag wires of dimensions L = 780 fim, 
w = 0.1 — 0.2 /im and t = 47 nm on oxidized Si sub- 
strates using electron-beam lithography, thermal evap- 
oration from a 99.9999%-purity Ag source, followed by 
lift-off of the bilayer resist. All the wires studied were 
evaporated simultaneously to ensure similar microstruc- 
ture and resistivity. The 15 samples were divided into 
four batches - one that was kept pure, and three that 
were implanted with 2, 6, and 10 ppm of Fe impurities, 
respectively. SRIM simulations of the implantation at 70 
keV indicate that over 90% of the implanted ions stay 
in the Ag wires, with the rest going into the substrate 
[lfij . As none of the 6 ppm samples survived, we re- 
port data taken on a pure sample, and samples with 2 
and 10 ppm Fe impurities. Table I lists the sample pa- 
rameters. All samples were measured immersed in the 
mixing chamber of a dilution refrigerator with filtered 
leads. Four-probe resistance measurements were made 
using a lock-in amplifier with a ratio transformer to im- 
prove sensitivity |lfj . The voltage drop across the sample 
was limited to eV s < fc^T to avoid heating the electrons. 

Raw magnetoresistance (MR) data for the three sam- 
ples, at T = 1.8 K, are shown in the inset to Fig. 1. MR 
data for each sample were fit using the following pro- 
cedure. Because Ag has moderate spin-orbit scattering 
(t so k 40 ps for all of our samples), the MR is positive 
at low temperatures, when 3> t so . At higher tem- 
peratures, when < t so , the MR starts out positive 
but then turns around at a field scale B ss 20 mT. Data 
at several temperatures in this higher temperature range 
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0.76 14.6 
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780 


47 110 


3890 
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0.86 9.1 
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AgFe2 


780 


47 185 


2330 


146 


0.72 4.1 
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TABLE I: Geometrical and electrical characteristics of the 
measured samples at 1.2 K. L, t, and w are the sample di- 
mensions and R is the resistance. D is the diffusion con- 
stant determined from the resistivity and the Einstein rela- 
tion p _1 = e 2 DvF, with the density of states in Ag vf = 



1.03 x 10 J~ m . L ao is the spin orbit length extracted 
from the fits of the magnetoresistance to weak localization 
theory. L™ ax is the maximum coherence length measured at 
40 mK. Ci m p is the implanted Fe concentration. 



(which was different for each sample) were first fit with 
three free parameters: L,p = (Dr^) 1 / 2 , where D is the dif- 
fusion constant, L so = (Dtso) 1 / 2 , and the sample width, 
w. For each sample, these fits gave consistent values of 
L so and w over a broad temperature range |lfj| . We then 
fixed those values of L so and w, and fit the MR curves 
for all temperatures with as the only free parameter. 
Using the value of D obtained from the resistance and 
sample dimensions, we finally obtain as a function of 
temperature for each sample. 
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FIG. 1: (color online) Total decoherence rate vs. temperature 
for a pure Ag sample (A), and for samples with implanted Fe 
concentrations of 2 ppm (•) and 10 ppm (■). Inset: Raw 
magnetoresistance data at T=1.8 K for the three samples, 
and fits to weak localization theory. 



Fig. 1 shows the decoherence rate, r^" 1 , obtained from 
the MR data and fitting procedure described above, for 
the three samples. The r^ 1 data for the pure sample 
follow theoretical expectations above 200 mK, with a 
modest amount of saturation at lower temperature |lfj| . 
Above 200 mK, the data are fit well (solid line) by a 
combination of electron-phonon scattering {t^ 1 oc T 3 ) 
and electron-electron scattering (t^ 1 cx T 2 / 3 ). The pref- 
actor of the T 2 / 3 term is equal to 0.36 ns _1 if~~ 2 / 3 , in 
good agreement with the theoretical value ^t|- The ad- 
ditional contribution to t^ 1 from the Fe impurities, y m , 
is clearly observable in the 2 ppm sample up to 4 K and 
in the 10 ppm sample up to 10 K. Above 10 K extraction 
of 7 m from t^ 1 is not reliable, due to an apparent differ- 
ence in the magnitude of the electron-phonon scattering 
rate in the 10 ppm and pure samples. 

The suppression of the weak localization magnetoresis- 
tance by magnetic impurities was first discussed in L9] 
for the case of static impurity spins. Later, Fal'ko lfj] 
pointed out that the impurity spins have their own in- 
ternal dynamics given by the Korringa rate, 7k °c ksT. 
Since spin-flip scattering of conduction electrons and re- 
orientation of the impurity spins arise from the same 
processes, the ratio of their rates is given by j m /jK = 
n s /upkBT, where n s is the density of magnetic impuri- 
ties and vp is the total density of states of conduction 
electrons at the Fermi level. Below Tk the Korringa rate 
saturates, so the ratio becomes "f m /"fK — tisIvf^bTk- 
If "fx > 7m j then the impurity spin configuration is 
randomized during the time between spin-flip scattering 
events of the conduction electrons. The contribution of 
magnetic scattering to the total decoherence rate in the 
spin-singlet channel of the weak localization magnetore- 
sistance is given by: 



27r, 



for 7 m > 7k. (1) 



for "f m <lK- (2) 



where is the dephasing rate due to electron-electron 
and electron-phonon scattering. For our Ag sample with 
2 ppm Fe, assuming Tk in the range 2-4 K (see below), 
we find "fm/lK = 0.02 — 0.04 <C 1, hence we should use 
Eq. J2J to extract 7 m from t^T 1 . Fig. 2 shows 7 m for 
the 2 ppm and 10 ppm samples obtained from the t^ 1 
data in Fig. 1 using Eq. J2J), with r-^ 1 obtained from the 
data on the nominally pure sample. The data sets from 
the 2 ppm and 10 ppm samples are consistent with each 
other when we use Eq. (J2Jl to analyze both of them. 

The criterion j m /jK < 1 is a necessarycondition for 
the theoretical approach of Micklitz et al. Jl2( in their nu- 
merical renormalization group (NRG) calculations of r y m . 
Those authors note, however, that Eq. is not strictly 
correct, because the decoherence induced by electron- 
electron interactions (l7] | is not an exponential process. 
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FIG. 2: (color online) Inelastic scattering rate due to magnetic 
impurities for the 2 ppm (•) and 10 ppm (■) samples. Data 
for the 2 ppm sample are multiplied by 5. The solid line is 
the theoretical calculation of Micklitz et al. fit to the 

data for T > 0.4 K. The dotted line is the Suhl-Nagaoka 
approximation for s=l/2. 



the theoretical estimate. (Those data were analyzed with 
the Suhl-Nagaoka approximation, which is in close agree- 
ment with the theory of for T > T K .) It thus appears 
that the inelastic scattering cross section of Fe in Ag is 
roughly twice that of Mn in Ag. 
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Eq. (12) in |l2j shows how to add the decoherence rate 
due to e-e interactions with the rate due to all other pro- 
cesses including j m . We have compared the results of 
analyzing our data with Eq. J2J and with Eq. (12) in 
[l2^ . and we find that for the 10 ppm sample the differ- 
ence is negligible. For the 2 ppm sample, using Eq. (12) 
in ^3] would increase 7 m by at most 15% at the lowest 
temperature, which is not significant on the logarithmic 
plot in Fig. 2. The solid line in Fig. 2 shows a fit of the 
NRG calculation from to our data at temperatures 
above 0.4 K only. Fitting to the 2 ppm (10 ppm) data 
alone gives Tk — 4.8 K (5.4 K). For comparison, the 
Suhl-Nagaoka approximation with Tk = 4.8 K is shown 
as the dotted line. The data and fits shown in Fig. 2 
lead immediately to several striking conclusions. First, 
as noted recently [l3| . the Suhl-Nagaoka approximation 
does not come close to reproducing the temperature de- 
pendence of 7 m for T < Tk- Second, the NRG theory of 
Micklitz et al. ,12} fits the data reasonably well over the 
temperature range T/T K = 0.1-2. And third, the NRG 
theory deviates strongly from the data when T/Tk < 0.1. 

According to ^2], the maximum value of 7 m , occur- 
ring at T = T K , is 0.23 x 4n s /(nhis F ). For 10 ppm 
of magnetic impurities in Ag, the theoretical estimate 
is T™ ax — 16 ns _1 , whereas the data in Fig. 2 show a 
maximum value about twice as large. This significant 
discrepan cy m ay point to the inadequacy of the spin-1/2 
theory of |l2j to account for the large spin (s=2) of Fe 
impurities in Ag. On the other hand, previous measure- 
ments of 7 m in Ag samples implanted with Mn impurities 
[To) , with s=5/2 and Tk ~ 40 mK, were consistent with 



FIG. 3: (color online) Change in resistivity vs. magnetic field 
for the 10 ppm sample at temperatures of 83, 188, 550, 860, 
and 1880 mK, from top to bottom. The zero for the 83 mK 
data is chosen so that the fit curve approaches Ap — » for 
\B\ 3> Bk- Subsequent data sets are shifted downward by 
0.75 nS7cm for clarity. Solid lines represent the global fit to 
the spin-1/2 NRG theory from |2cj|. as discussed in the text. 
Left inset: Raw Ap(B) data at 83 mK and 10 K. The latter 
shows the classical magnetoresistance Ap{B) oc B 2 , which 
was subtracted from the raw data sets to produce the curves 
in the main figure. Right inset: Ap(B = 0) vs. temperature. 

To further test the theory, we seek an independent es- 
timate of Tk to compare with the value obtained from 
the fit to 7 m shown in Fig. 2. The traditional method 
of determining Tk from the resistivity vs. temperature, 
p(T), is not reliable for our samples, because the electron- 
phonon scattering contribution to p(T) starts to grow sig- 
nificantly when T > 8 K, before the Kondo contribution 
has completely died out. Instead, we analyze the high- 
field magnetoresistivity, p(B), shown at several temper- 
atures in Fig. 3. p(B) data for T = 53 mK (not shown) 
are indistinguishable from those at 83 mK, indicating 
that the Kondo contribution to p is close to the unitar- 
ity limit at those temperatures. We fit the p{B) data 
using NRG calculations from ref. Hj. For T <C T K , 
these are nearly identical to the T = analytical solu- 
tion from [2l| . The two-parameter fit to the 83 mK data 
set gives the values ApK = p{B = 0, T = 0) — p(B 3> 
Bk, T = 0) = 0.57 nOcm/ppm and B K = 1.2 T, where 
gpsBK = IcbTk, with ps the Bohr magneton |22|. Us- 
ing g=2, that in turn gives Tk = 1.6 K. Similar fits to 
the p(B) data sets at higher temperatures, keeping ApK 
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fixed, show that they can not be fit consistently with such 
a low value of Tk ■ 

To circumvent this problem, we have performed a 
global fit to all the p(B) data sets with only three pa- 
rameters: ApK, Bk, and Tk, where Bk determines the 
magnetic field scale over which p(B) decreases, while Tk 
determines the temperature scale over which p{B = 0) 
decreases. The results of this global fit are shown by 
the solid lines in Fig. 3, with the values Bk = 1-36 T, 
T K = 2.96 K, and Ap K = 0.58 nOcm/ppm. The value of 
Bk corresponds to a Kondo temperature of 1.8 K. The 
discrepancy between the values of Tk extracted from the 
field scale and from the temperature scale highlights the 
fact that we are using a spin-1/2 theory. For free spins 
in a metal without Kondo effect, the magnetoresistance 
is proportional to — < s z > 2 , where < s z > is the aver- 
age spin polarization given by the appropriate Brillouin 
function. A comparison of the functions — < s z > 2 for 
s=l/2 and s=2 shows that the full width at half max- 
imum of the s=l/2 function is 1.85 times larger than 
for the s=2 function. If a similar relationship holds for 
the p(B) curves in the Kondo regime, then the value 
B K = 1.2T we found from the s=l/2 fit to the 83 mK 
data would become Bk ~ 2.2 T for s=2, which corre- 
sponds to a Kondo temperature Tk = 3.0 K. That is 
consistent with the value of Tk obtained from the tem- 
perature dependence of p(B = 0), but still lower than 
the value Tk ~ 5 K obtained from the fit to j m . 

The experimental value ApK — 0.58 nf2cm/ppm 
should be compared with the unitarity limit for s-wave 
scattering, ApK = 4:irHn s /ne 2 kF — 0.43 nficm/ppm 
for Ag, where Uf is the Fermi wavevector. (The value 
for d-wave scattering is five times higher.) Our mea- 
sured value of ApK is larger than the theoretical s-wave 
value, but less than the measured values in other Fe- 
noble metal Kondo systems: 1.3 nflcm/ppm in Cu:Fe 
and 1.0 nf2cm/ppm in Au:Fe [23j |. 

The most serious discrepancy between theory and ex- 
periment is the flattening of 7 m (T) for T/Tk < 0.1, 
shown in Fig. 2. The behavior of the data is reminiscent 
of Eq. (3b) in [24|], which describes inelastic scattering by 
an underscreened Kondo impurity when T <C Tk- Since 
the magnitude of j m at T — 40 mK is 40 times smaller 
than at 4 K, however, fitting the low-T data with Eq. 
(3b) of [24| would require a substantial reduction of the 
prefactor. The issue of screening for Fe in Ag is unclear. 
At first glance, the five channels corresponding to the 
d-orbitals should overscreen the s=2 Fe spin. In reality, 
however, the different channels have different coupling 
constants, and the impurity spin and orbital degenera- 
cies are broken by crystal fields and spin-orbit scattering. 
Further theoretical work on realistic models of spin-2 sys- 
tems should help to resolve these issues. 

Note added: We recently became aware of a related 
work leading to similar conclusions |25| . 
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